Pionic decay of A Hypernuclei in a Continuum Shell 

Model 



C. Albertus, J.E. Amaro and J. Nieves 

Departamento de Fi'sica Moderna, Universidad de Granada, E-18071 Granada, Spain. 



Abstract 

We evaluate pionic decay widths of A hypernuclei using a shell model for, both 
the nuclear bound and the continuum nucleon wave functions in the final state, and 
distorted waves for the outgoing pion. An excellent agreement with the recent KEK 
measurement of 7r~-decay widths of and 2 ^Si is found. Besides, results for 5 ^Fe 
are consistent with the existing upper bound. 
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1 Introduction 



The pionic decay of A hypernuclei has received attention in the past both from the exper- 
imental [l]-[7] and theoretical [8]-[22] points of view (for a recent review see Ref. [23]). 
It is well known that the mesonic decay is largely suppressed by Pauli blocking, although 
the consideration of proper pion distorted waves [13, 16, 20] weakens the effect of the 
Pauli blocking considerably. On the other hand, it was also pointed out in Ref. [16] the 
relevance of the use of a correct energy balance in the decay. Two theoretical frameworks 
have been traditionally used to describe this process: Polarization Propagator Method 
(PPM) [12, 18, 21], supplemented by quantum field theory functional techniques [22, 23], 
and the Wave Function Method (WFM) [13]-[17]. In the PPM, the calculation of the 
pionic widths is performed in nuclear matter and the results are translated to finite nuclei 
by means of the Local Density Approximation (LDA). Both methods PPM and WFM 
turn out to be related in a semiclassical limit [19], being the WFM more reliable than the 
PPM+LDA, since the mesonic decay channel is quite sensitive to the details of the shell 
structure of the hypernuclei, specially for light systems. 

In Ref. [16] a simple Woods-Saxon (WS) shell model was employed, with a central 
potential of constant depth for several nuclei. This global potential did not include spin- 
orbit term and only the radius of the WS well depended on the specific nuclei (oc A3). 
The resulting binding energies of the shells were globally shifted, by hand, to reproduce 
the ground state masses of the involved nuclei. Besides, the continuum contribution to 
the decay was estimated by discretizing the positive nucleon energy levels by means of an 
infinite barrier placed at distances of about 20 fm. Despite all these crude approximations, 
the model led to predictions for the mesonic decay widths of 12 A C which were compatible, 
within errors, with the available experimental measurements at that time, for both 7r° and 
7r~ decay channels. Furthermore, the model of Ref. [16] provided an overall description 
of medium and heavy hypernuclei mesonic decay, which is still nowadays accepted as one 
of the most reliable theoretical estimates [23] . 

Very recently, precise measurements of the ir~ mesonic decay of 12 A C and 2 ^Si have been 
obtained at KEK, and an upper bound for 5 ^Fe has also been given [7]. The purpose of 
this paper is to update the model of Ref. [16], improving the nuclear structure description 
of the process, by using: i) more realistic potentials for each one of the nuclei by fitting 
the parameters of a WS well, with spin-orbit forces, to the ground state and also to some 
of the discrete final state energies of the nuclei involved in the decay and ii) a Continuum 
Shell Model (CSM) to describe the positive energy tail of the mesonic decay width. In 
this way we replace the discrete sum used in Ref. [16] by the appropriate integral over the 
continuum states. This correct treatment of the continuum contribution is important, 
since, as we shall see, for 7r~- 2 ^Si decay it amounts about 50% of the total measured 
mesonic width, in contrast to the case of the hypernuclei studied in Refs. [13] and [16], 
for which the continuum contribution was much smaller. 

Since, as mentioned above, the use of pion distorted waves turns out to be crucial [16], 
[20], the improved treatment of the nuclear structure allows us to use the new and precise 
measurements performed at KEK, as a further test of the quality of the pion-nucleus 
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dynamics used in Ref. [16] and developed in the works of Refs. [24] and [25]. A deep and 
detailed knowledge of the pion dynamics inside of a nuclear medium has become a topic 
of renewed interest to explore possible partial chiral restoration in the medium [26]- [28]. 

The paper is organized as follows: In Sect. 2 the needed formulae to compute the 
pionic decay width are given, both for the discrete and continuum contributions. Details 
on the nuclear CSM used in this work are discussed in Sect. 3, where the energy reaction 
balance is also studied. Finally in Sect. 4 we present the results of this work and our main 
conclusions. 



2 Theoretical Description of the Pionic Decay. 

In this work, we compute the mesonic decay width as the sum of the contributions of the 
following processes 1 

z -+ (*z) d + A (i) 

Z - (^Z + I^ + tt- (2) 
Z - +n + vr°, (3) 

V / gs 



where d denotes the ground or discrete excited states of the final nucleus. As we shall see, 
we evaluate the processes of Eqs.(l), (2) by putting the outgoing nucleon, coming from 
the weak A decay, in an unoccupied bound shell of the daughter nucleus. On the other 
hand, in Eqs. (3) and (4) gs means that the daughter nucleus is left on its ground state. 
In the two last reactions the outgoing nucleon, coming from the weak A decay, goes to the 
continuum (positive energy) and we denote that contribution by T c , while the first two 
reactions give what we call discrete contribution T d . Thus we split the pionic width into 
two contributions, T = T c + IV Experimentally, what can be observed are the inclusive 
processes 



A 



Z -> X + 7T°, (5) 

Z -> X + 7i-. (6) 



The main contribution to these processes is given by the exclusive ones shown in Eqs. (1- 
4). 

The pionic decay is produced by a baryonic one-body operator 

5Hl N = -Gml \S-— a- q cm } r\ (7) 

where (Gml) 2 /8n = 1.945 x 1(T 15 , the constants S and P are equal to 1.06 and 0.527 
respectively and is the pion mass (139.57 or 134.98 MeV for 7r~ or 7r°), q cm is the 

1 We use the notation A Z, to denote a nuclear core A ~ 1 Z and one bound A particle. 
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momentum of the outgoing pion in the A rest frame, and the Pauli matrices a and r A , 
where A is a cartesian isospin index which will be contracted to the pion field, act on the 
spin and isospin Hilbert spaces respectively. Taking the A isospin wave function as that 
of a neutron, the r A operator in Eq. (7) implements the extreme AT = 1/2 rule, which 
leads to a rate of A — > ix~p twice as large as that of A — > n°n. 

The free space A decay width is readily evaluated and leads for proton or neutron 
decay to 

q cm 2Ma , w 

X(x,y,z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz, (10) 

where a indicates neutron or proton or equivalently it or 7T~ decay, C ^ =4 and C ^ = 2 
are isospin coefficients, and M (938.27 or 939.57 MeV for p or n) and M A (1115.68 MeV) 
are the nucleon and A masses respectively. The total free space hyperon A width, I\, is 
given by the sum of proton (ir~) and neutron (ir°) contributions, 

r A = + rf"] e (ii) 



In the case of a bound hypernucleus, assuming that the hyperon A is in the lsi shell 
and a closed shell structure for the underlying nuclear system, A_1 Z, the width for any of 
the processes of Eqs. (1-4) is given by 2 

x (Gml) 2 

where nljm and stands for the quantum numbers and relativistic energy of the outgo- 
ing nucleon in the nlj shell, |A, m SA ) denotes the A wave function with angular momentum 
third component m SA , the A energy, including its mass, u(q) the pion energy, and the 
sum over N runs over the unoccupied nuclear orbitals (n, I, j, m). In Eq. (12) the 
sums over N are over proton or neutron orbitals according to a. The pion wave function 
(g, x)*) as a block corresponds to an incoming solution of the Klein Gordon equation, 

^\q,x)* = (uj(q) - V c (x)) 2 ffi(q,x)*, (13) 



A, m 



SA 



S a- 



nljm;ip^\q) 



(12) 



-V + ml + 2u(q)V opt (x) 



with Vc(x) the Coulomb potential created by the nucleus considering finite size and 
vacuum polarization effects, for tt~ and zero for n°, and V opt (x) the optical potential which 

2 The following expressions are valid for the discrete contribution (Eqs. (1) and (2)) to the decay width. 
The needed modifications to compute the continuum part (Eqs. (3) and (4)) are discussed after Eq. (24). 
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describes the 7r-nucleus interaction. This potential has been developed microscopically 
and it is exposed in detail in Refs. [24, 25]. It contains the ordinary lowest order optical 
potential pieces constructed from the s- and p-wave nN amplitudes. In addition second 
order terms in both s- and p-waves, responsible for pion absorption, are also considered. 
Standard corrections, as second-order Pauli re-scattering term, ATT term, Lorentz-Lorenz 
effect and long and short range nuclear correlations, are also taken into account. This 
theoretical potential reproduces fairly well the data of pionic atoms (binding energies and 
strong absorption widths) [24] and low energy 7r-nucleus scattering [25] . 

After a little Racah algebra, one gets for the A decay width inside of a nucleus 



r («) = £ r H 

N=nlj>F 

-p(a) _ ^ 
1 N — 



An 



\Gmif- 



Qn 



+ u(q N )/E A 



with 



q N = ((E A - E. 

7 « 



N 



2x1/2 



(14) 

(15) 
(16) 



Note that E^,qN and the integrals S^' and defined below, depend on the isospin 
a. We have implemented the recoil factor (1 + ujj ' Ea)~ 1 , being Ea the energy (including 
the mass) of the daughter nucleus, because most of the decay corresponds to nucleons in 
nuclear bound excited states, and as a consequence the nucleus of mass Ma recoils as a 
whole. S$(qN) and S^^n) are the s-wave and p-wave contributions given by 



S nlj{QN) 



2j + 1 , r , , |2 



N 



4a U WniMN)? hi-1/2 + + l)\Nm,(q N )\ 2 S jtl+1/2 } 
Qn 



(17) 
(18) 



with 



Inlj(qN) 

M nlj (q N ) 
N„ij{q N ) 



r>drR[ k Xr)Rr{qN,r)R r a 3 {r) 



f rHrRt\r) ( ^^1 - (/ _ 

fVdr<>(r) ( d Jm^l + (Z + 2) <W) 
Jo \ dr r 



(19) 

Rnij(r) (20) 
Rm 3 {r). (21) 



Here R\ w \qN',r) are the radial wave functions of the pion for each partial wave, regular 
in the origin and with the asymptotic behavior 



1 7T 

~ e lSl — sin(gr — I— + 5i) for n° 
qr ' 2 



(22) 
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Rl n \q; r^oo ~ e i{Si+ai) — sin(gr - I- + o x + <J, - 77 log2gr) for tT (23) 

t^T* 2 

with 77 and 07 the Coulomb parameter and phase shift defined as in Ref. [29] and Si the 
complex (to take into account inelasticities) phase shifts obtained from the numerical 
solution of the Klein Gordon equation. Finally, and R n ij(r) are the A and nucleon 
bound radial wave functions, normalized as usual 

r+00 

/ drr 2 \R(r)\ 2 = 1 (24) 
Jo 

The A wave function in the initial hypernucleus, R^ A \ is obtained from a WS potential [33] 
to account for the mean A-nuclear core interaction, with parameters compiled in Table 1. 
The nucleon dynamics will be studied in detail in Sect. 3. 

Some of the nuclei which we use are not closed shell nuclei. In this case the nucleons 
from the A decay can fill up empty states in a n, l,j shell. We take that into account 
by multiplying S$ and S$ by n h /(2j + 1). 



Hypernucleus V A [MeV] R A [fm] a A [fm] B A [MeV] 



12 p 


31.1 


2.45 


0.60 


-10.8 [30] 


2 isi 


30.1 


3.30 


0.60 


-16.6 [31] 


56 A Fe 


30.3 


4.18 


0.60 


-21.0 [11] 



Table 1: Lambda hyperon WS parameters and ground state (lsi) binding energies (B A ). The total A 
energy, E A is given by M A + B A . 



The above equations (12)-(24) can be readily used to compute the discrete contribution 
to the pionic decay width, but when the outgoing nucleon coming from the weak A decay 
goes to the continuum (positive energy), one should replace the sum Y,N=nij>F m Eq. (14) 
by a sum over multipoles, plus an integral over the nucleon energies, i.e., 

Rnij(r) -> Rij(p;r) (26) 



with E and p = y/ E 2 — M 2 the nucleon energy and momentum, E ma _ x = E A — m n , the 
maximum nucleon energy, neglecting the recoil of the nucleus, and Rij(p;r) a continuum 
solution of the nucleon Schrodinger equation [32], with the same potential as that used 
for the bound nucleons to compute the discrete contribution to the decay width. The 
nucleon wave function normalization in the continuum is the same as that given above 
for the pions (Eqs. (22-23)). 

As a test of the multipole expansion in the continuum, we recovered the free space A- 
decay width of Eq. (8) from the multipolar expansion of Eqs. (14) to (26). For this purpose 
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1 and replaced the radial pion and nucleon wave functions by spherical Bessel's 

47rL 3 /3 the volume 



we set A 

functions, ji, and the radial A-wave function by \jAn/V, being V 
of interaction, which will be finally sent to infinity when calculating physical observables. 
The test is straightforward, taking into account: 



1=0 



i — 2 

/47T f L 

V ~v lo drr • 7i ( pr )^ r ) 

f^drr 2 {^{21 + l)jf(pr)^ £ drr 2 3l (pr) 3l (qr) + 0(1/ L) 



7T 



= ^(P - ?) + 



(27) 



On the other hand, we have also used this free space limit to test our computational code 
and the precision of our numerical algorithms to solve differential equations, to perform 
integrations and the sum over an infinite number of multipoles. The idea was, taking as 
starting point the formulae of Eqs. (14) to (26) with a radial A- wave function given by 



An/V and switching off the pion-nucleus, and the nucleon potentials, to recover numer- 
ically, in the limit L going to infinity, the free space A decay width. Such a calculation, 
from the numerical point of view, is much more demanding that the actual calculation of 
the pionic decay width of the hypernuclei, since the Dirac's delta appearing in Eq. (27) has 
to be constructed numerically out of a large sum over multipoles and slowly convergent 
integrals. 



3 Nuclear Structure and Energy Balance 

We model the nuclear structure of the Z system by a Slater determinant built with 
single-particle wave functions obtained by diagonalizing a WS potential well 

V ws (r) = V f(r, R , a ) + V J^^Il^A + V c (r), (28) 

r dr 



where 



f(r,R,a) = - (29) 

1 + exp 



and Vc(r) is the Coulomb potential created by an homogeneous charge distribution of 
radius Rc- The parameters, compiled in Table 2 for neutrons and Table 3 for protons, 
of this potential are adjusted 3 to reproduce some experimental single-particle energies 
around the Fermi level [34]- [35]. This will enforce not only a correct energy balance 
for the decay process to the first available shell but also to some excited states. The 
main contributions to the processes of Eqs. (1) and (2) come from situations where the 

3 Essentially, we adjust the depths V and V^s, and for the radius- and thickness-type parameters 
standard values have been used 
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V [MeV] 


Ro H 


a [fm] 


V LS [MeV] 


Rls [fm] 


a LS [fm] 


11 C 


-57.0 


2.86 


0.53 


-6.05 


2.86 


0.53 


27 Si 


-66.2 


3.50 


0.70 


-3.30 


3.75 


0.70 


55p e 


-54.0 


4.70 


0.50 


-8.30 


4.70 


0.50 



Table 2: Neutron WS parameters used in this work. 





V [MeV] 


Bo [fm] 


a [fm] 


V LS [MeV] 


Rls [fm] 


a LS [fm] 


Rc [fm] 


11 C 


-38.4 


2.86 


0.53 


-6.05 


2.86 


0.53 


2.86 


27 Si 


-47.4 


3.75 


0.53 


-10.0 


3.75 


0.53 


3.75 


55p e 


-50.8 


4.70 


0.50 


-8.30 


4.70 


0.50 


4.70 



Table 3: Proton WS parameters used in this work. 

daughter nucleus is left in the ground state or in the first few excited states. Since the 
effect of the Pauli blocking depends strongly on the pion energy after the decay it is 
important, as shown in Ref. [16], to perform a correct balance of energies, using when 
possible experimental energies. 

Thus, the energy of the first non-occupied shell is fixed to the mass difference between 
the ground states of the A Z and A ~ X Z nuclei for the case of neutron (ir°) decay and of 
the A (Z+ 1) and A ~ X Z for the case of the proton (rr ) decay channel. In Table 4 we give 
the first available shells and their empirical energies obtained in this way. The energies 
of this table do not totally determine the parameters of the mean field nucleon potential, 
and one has still the possibility to fit some excited state energies. In what follows, we 
give some more details of the adjusted shells. 



Hypernucleus 


Neutron Shell Energy [MeV] 


Proton Shell Energy [MeV] 


12 p 

2 isi 

5 A Fe 


lp 3 /2 -18.72 
ld 5 /2 -17.18 
2p 3/2 -11.20 


lpi /2 -0.60 
2si /2 -2.07 
I/7/2 -5.85 



Table 4: Single particle energies for the first available, to the pionic decay of the hypernuclei studied 
in this work, nucleon shells and their empirical binding energies. Energies have been obtained from the 
neutron and proton separation energies of the nuclear A Z and A (Z + 1) spices respectively, and have been 
taken from Ref. [36]. 
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• \C: We assume a closed proton and neutron 1^3/2 shell configurations for the 
ground state of 12 C and a lpi/2 (proton) lp^} 2 (neutron) configuration for the ground 
state of 12 N. The rest of neutron and proton shells adjusted by the mean field 
potentials, which parameters are given in Tables 2 and 3, are: 

Neutrons: We fix the energy of the lpi/2 neutron shell from the energy of the first 
excited state of 12 C [2 + 12 C*(4.43) MeV]. We assume a particle-hole configuration 
lpi/2, IP3/2 which leads to a binding energy for the lpi/2 shell of —14.29 MeV. 
We are aware that this excited state of 12 C might also have some contributions 
from proton degrees of freedom, configuration mixing, 2p2h or higher excitations, 
collective degrees of freedom, etc, not considered in this simple model. However, as 
long as this state had a sizeable component | lpi/2, 1^3/2)5 ^ would be reachable in the 
weak decay of the hypernucleus and the procedure described above would guarantee 
that the energy balance is correct, not only when the 12 C is left on its ground state, 
but also when it is left on its first excited state, leading then in both cases to good 
pion wave-functions. The nuclear matrix elements appearing in the evaluation of 
the decay width, though depending strongly on the pion wave-function, are less 
sensitive to the specific details of the nuclear wave function. 

Protons: There are no 12 N excited states amenable to be explained within a shell 
model, since all levels compiled in Ref. [36] are broad resonant states populated in 
nuclear reactions. Therefore, we choose to use the same spin-orbit force as in the 
neutron case, and fix the depth of the central part of the potential to reproduce the 
lpi/2 shell as given in Table 4. The mean field potential adjusted in this way, does 
not provide excited states for 12 N. 

• 28 \Si: We assume closed proton and neutron ld^/2 shell configurations for the ground 
state of 28 Si and a 2s±/2 (proton) ld^} 2 (neutron) configuration for the ground state 
of 28 P. The rest of neutron and proton shells adjusted by the mean field potentials 
are: 

Neutrons: We fix the energy of the 2si/ 2 neutron shell from the energy of the 
first excited state of 28 Si [2 + 28 Si*(1.78) MeV]. We assume a neutron particle-hole 
configuration of the type \2si/ 2 , lrf^A) which leads to a binding energy for the 2s!/ 2 
shell of 15.40 MeV. Limitations and virtues of describing the mesonic decay with 
this simple picture for the underlying nuclear core, are similar to those commented 
above for the l \C case. We should mention that to adjust the WS neutron potential 
to simultaneously give the empirical I0I5/2 and 2si/ 2 energy shells, keeping the ld^/2 
energy shell above the 2si/ 2 , is delicate and that we had to use a value a bit high 
(0.70 fm) for the thickness parameters Oq and a^s- We have also tried smaller values 
for the thickness. For instance for values of a = Ols = 0.58 fm, in order to adjust 
the empirical energies of both the ld 5 / 2 and 2si/ 2 shells, a spin orbit force very small 
(y LS = -0.1 MeV with V = -65.75 MeV) is required. Despite that both the ld 5/2 
and ld 3 /2 shells turn out to be almost degenerate, and that the latter one is now 
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deeper than the 2s\/ 2 shell, the decay width is rather stable, and it gets increased 
only by about 10%, with respect to the results of the next section. This increase 
is due to an enhancement of the shell contribution, but the corresponding 

configuration \ld,3/ 2 , 1^5/2) has an excitation energy too small to correspond to any 
experimental excited state of 28 Si. 

Protons: For the 28 P nucleus, as for the case of 12 N, there are no excited states 
amenable to be explained within a shell model, since all levels compiled in Ref. [36] 
are again broad resonant states populated in nuclear reactions. We fix the depth 
of the central part of the potential to reproduce the 2s\j 2 shell as given in Table 4 
and we choose the depth of the spin-orbit such that the next shell (1^3/2) appears 
in the continuum. We will discuss in the next section, the dependence of our results 
on the precise value of the spin-orbit force. If we use instead the same spin-orbit 
force as in the neutron case, there will appear excited 28 P states, which cannot be 
identified in the experiment. 

• 56 A Fe: The nuclear core structure of this hypernucleus is more difficult to describe, 
within our simple shell model, than those of 28 v Si and 12 A C hypernuclei. Thus, our re- 
sults for the decay of this hypernucleus are subject to more theoretical uncertainties. 
In what the ground states respect, for 56 Fe, we assume a configuration composed of 
two paired I/7/2 proton holes and two paired 2j>3/2 neutron particles, while for 56 Co, 
we assume a proton hole in the I/7/2 shell and a neutron in the 2p 3 / 2 shell. 

The rest of neutron and proton shells adjusted by the mean field potentials are: 

Neutrons: We fix the energy of the I/5/2 neutron shell from the energy of the first 
excited state of 56 Fe [2 + 56 Fe* (0.847) MeV]. We assume for neutrons a two particle 
configuration |2p 3 / 2 , 1/5/2) which leads to a binding energy for the I/5/2 shell of 
10.35 MeV. 

Protons: The first excited state of 56 Co has spin-parity 3 + and an excitation energy 
of 0.16 MeV. In principle, one might use it to determine properties of the WS proton 
mean field potential. However, a word of caution must be said here. One might try 
to describe this state as a proton configuration with two paired holes in the 1/7/2 
shell and a particle in the 2p 3 / 2 shell, in such a way that the above configuration 
would determine the energy of the 2p 3 / 2 shell. The I/7/2 shell completes 28 protons 
which is a magic number, and thus one expects an energy gap between this shell 
and the following in energy, 2p 3 / 2 , appreciable and of the order of the MeV, and not 
as small as 0.16 MeV. Then, likely, the 3 + excited state should have an important 
neutron component (one neutron in the 2p 3 / 2 shell and the other in the I/5/2 one) or 
more complex components not considered in our simple shell model. Then, it seems 
safe to guarantee that this state will not have a sizeable component \2p 3 / 2 ,lfy^) . 
Therefore, we choose to use the same spin-orbit force as in the neutron case and fix 
the depth of the central part of the potential to reproduce the 1/7/2 shell as given 
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7T° Decay 


7r Decay 


V op t 


r c r r cxp [3] 


r d r c r r exp [7] 


12 n FP 

aL NQ 


0.136 0.008 0.144 

0.217 ± 0.084 

0.150 0.008 0.158 


0.079 0.027 0.106 

0.113 ±0.013 ±0.005 

0.082 0.028 0.110 


FP 

2 A Si 

A NQ 


0.061 0.001 0.062 
0.074 0.001 0.075 


0.018 0.019 0.037 

0.047 ±0.008 ±0.002 

0.020 0.019 0.039 


A Fe 

A NQ 


0.010 0.003 0.013 
0.010 0.003 0.013 


0.005 0.009 0.014 , M 

< 0.015 90%C.L. 
0.004 0.010 0.014 v ; 



Table 5: Pionic decay widths, units of T\, calculated with two different pion nucleus optical potentials: 
FP stands for the full optical potential of Ref . [25] , NQ stands for the pion- nucleus interaction obtained 
by switching off the imaginary part of the FP optical potential coming from quasielastic pion scattering. 

in Table 4. The mean field potential adjusted in this way, leads to an excited state 
of about 5 MeV above the ground state. 

4 Results and Concluding Remarks 

Results for the pionic decay widths of 12 \C, 2 \Si and 5 \Fe, calculated with two different 
pion nucleus optical potentials are presented in Table 5. The effect of the imaginary 
part of the potential is to remove from the emerging pion flux those pions which undergo 
quasielastic scattering or pion absorption. However, while the pions absorbed should be 
definitely removed, this is not the case with those which undergo quasielastic scattering, 
since even if they collide, they are still there and will be observed. This means that one 
should not remove these pions from the pion flux and we take this into account here. The 
effect is moderately small, as it was already noted in Ref. [16]. 

The agreement with the recent KEK measurements is remarkably good and it is also 
good when our results are compared to the older measurement of the n° decay width 
of 12 A C. This is a clear success of the model of Refs. [24, 25] to account for the pion- 
nucleus dynamics at low energies. In Ref. [16], it was obtained a value of 0.086 T A for 
the 7r~- 12 A C decay width. In this work we find a value about 25% higher and in a closer 
agreement with the experiment. Differences are even bigger if one looks at the continuum 
contribution and also appear for the n° decay case. In both works, here and in that of 
Ref. [16], the same 7r~-wave function has been used, being then the difference due to an 
improved treatment of the underlying nuclear core dynamics. As we will see below, the 
barrier method employed in Ref. [16] to estimate the continuum contribution compare 
reasonably well to the more correct treatment followed here, when the same nuclear 
potential is used. The discrepancies have to be attributed not only to the different WS 
potentials used in both works, but also to the somewhat artificial procedure followed in 
Ref. [16] to enforce the correct energy balance in the decay. 
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E-M [MeV] E - M [MeV] 

Figure 1: Continuum proton energy distributions from the ir~- decay width of 2 ^Si. Results have been 
obtained with the NQ 7r~-nuclcus optical potential, as defined in the caption of Table 5, and only the 
most relevant multipoles are shown. 



Note that in Table 5 the computed widths are always slightly below the data. The 
contribution of other reaction channels, not included in Eqs. (1-4), though we expect it 
to be small, would improve the present calculation. 

In contrast to most of the hypernuclei studied previously by using the WFM [13]-[17], 
for the three hypernuclei considered here, the continuum contribution in some cases plays 
a crucial role. As a matter of example, for 7r~ decay of 2 ^Si, it turns out that T c is of 
the same size as T d (0.020 vs 0.019) and essential to understand the experimental datum. 
This reinforces the need of updating the calculation of Ref. [16], where the evaluation of 
the continuum contribution was a bit rough, since there was assumed that it was only a 
small fraction of the total. Indeed, in Ref. [16], the continuum contribution was estimated 
by discretizing the positive nucleon energy levels, by means of an infinite barrier placed 
at distances of about 20 fm. In Table 6 we compare the CSM of this work with the model 
of Ref. [16]. In both cases we use the same nucleon WS mean potential. Both methods 
agree remarkably well and much better of what one might expect a priory. Within the 
model of Ref. [16], the sizeable contribution of the continuum is due to the presence of 
a quasi-bound (1.27 MeV) state, ld^/2 and it has its counterpart in the size of the d%/2 
multipole in the CSM. Even more, the differential partial width dT^ 2 /dE presents a 
narrow peak (resonance), which gives most of the contribution to the integrated partial 
width, located around 1.27 MeV, as can be seen in Fig 1. Small changes in the WS 
proton mean potential can bind this shell, going then its important contribution to the 
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discrete part, T^, but remaining the total width + T c almost unchanged. For example, 
if one uses a spin-orbit force depth of —7 MeV instead of —10 MeV, the lofa/2 proton 
shell becomes bound (—0.1 MeV) and the total width is 0.042 I\ instead of the value of 
0.039 I\ quoted in Table 5. In the model of Ref. [16], the exact position of the barrier 
might influence energies, the number of shells, and contributions of each shell, but again 
the total contribution to T c remains rather constant, as long as the barrier is placed far 
enough. 

We have tested for the sensitivity of our results to the mean field nucleon WS pa- 
rameters. Thus, we have increased and decreased the spin-orbit depth Vls by 10% and 
re-adjusted the depth of the central part of the nucleon potential, Vo, to get the experi- 
mental ground states masses of the involved nuclei, it is to say V is modified to reproduce 
again the energies given in Table 4, with the new values of the spin-orbit force. Results 
turn out to be quite stable, changing at most by a 2%, except for the n~- decay of 5 \Fe, 
where the uncertainty can be as large as one part in fifteen. 

To finish, we show the continuum nucleon energy distributions for the first nucleon 
multipoles, which give the bulk of the total, (Figs. 1 to 4) and the contribution to 
T d of each of the unoccupied shells (Table 7) for the 7i~ and 7T° decay widths of 12 A C, 
of 28 A Si and of 5 ^Fe. All results have been obtained using the NQ 7r-nucleus optical 
potential. Resonances appear as distinctive features in the continuum contribution of 
some multipoles. 



12 



Neutron Decay 



Proton Decay 




> 1.6 x 10 



< 
f-< 



^3 

S u 4 x 10 
f-< 



> 



0.00075 
0.0006 



< 0.00045 f- 




Et) 0.0003 - 

^3 



t-, 



> 6 x 10- 



< 
f-< 



t-, 




> 

v 0.00075 



< 



t-, 

"X3 



0.0005 



0.00025 




15 20 



— 0.01 



> 



— 0.0001 



Eh 



"a 



Eh 

^3 




5 10 15 20 
E — M [MeV] 



i 

> 



-si 

Eh 

^3 



0.003 - 



0.002 



0.001 





^3/2 







5 10 15 
E-M [MeV] 



20 



Figure 3: Neutron (left) and proton (right) energy distribution for 1 ^C-decay. Solid (dashed) line 
corresponds to the j = I + 1/2 (j = I — 1/2) multipolc. The neutron integrated decay widths arc 
1.4 x 1(T 4 , 3.0 x 1(T 4 , 7.1 x 1(T 4 , 6.62 x 1(T 3 , and 4 x 1(T 5 , in units of T A , for the ls 1/2 , lp 1/2 , 
lps/2, 1^3/2 an d 1^5/2 multipoles respectively. While for protons, the contribution of the multipoles are 
4.33 x 1CT 3 , 5.79 x 1CT 3 , 2.16 x 1CT 3 , 4.30 x 10~ 3 , and 1.124 x 10~ 2 . 
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Figure 4: The same as Fig. 3 for 5 ^Fe-decay. The neutron integrated decay widths are 8.2 x 10~ 4 , 
2.3 x 1(T 4 , 3.2 x 1CT 4 , 1.56 x 1CT 3 , and 1 x 1CT 5 , in units of T A , for the ls 1/2 , lp 1/2 , lp 3 / 2 , l«?3/2 
and ld 5 / 2 multipoles respectively. While for protons, the contribution of the multipoles are 2.47 x 1CP 3 , 
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This work 
Multipole 


Barrier at 20 fm [16] 
Shell Energy [MeV] [r,JJ] 


S1/2 0.0007 


3si /2 2.45 0.0002 
4si /2 5.19 0.0003 
5si/ 2 9.21 0.0002 
Total s V 2 0.0007 


p 3 /2 0.0011 


3p 3 /2 4.59 0.0005 
4p 3/2 7.25 0.0004 
Total p 3 /2 0.0009 


Pi /2 0.0007 


3pi /2 5.02 0.0002 
4p 1/2 8.14 0.0002 
Total pi/2 0.0004 


4/2 0.0165 


14/2 1-27 0.0160 
24/2 3.22 0.0001 
34/2 6.31 0.0001 
Total pi/2 0.0162 


T c = 0.019 


T c = 0.018 



Table 6: Continuum contributions to the pionic decay width, units of Ta, calculated with the NQ 
7r~-nucleus optical potential, as defined in the caption of Table 5, and two different methods: the sum 
of multipolcs of the type defined in Eqs. (14) and (25) and the sum over positive energy discrete 
bound (by the effect of an infinite barrier placed at 20 fm) states [16]. Only contributions to T c larger 
than 5 x 1CP 4 Ta are shown. Results are for the 7r~-decay of the 2 A Si hypernucleus. 
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A 7 



7T° Decay 
Shell Energy [MeV 



nlj 



7r Decay 
Shell Energy [MeV 



nlj 



12 ri 



1P3/2 


-18.72 


lPl/2 


-14.29 


14/2 


-3.51 


2si/ 2 


-2.07 




Total 


14/2 


-17.18 


2si/ 2 


-15.40 


14/2 


-14.76 


2^3/2 


-3.29 


1/7/2 


-2.95 


2pi/ 2 


-2.35 




Total 



0.0473 
0.0865 
0.0156 
0.0008 
0.150 



IP1/2 



-0.60 



Total 



0.0817 



0.082 



28 A Si 



0.0086 
0.0285 
0.0342 
0.0009 
0.0010 
0.0009 
0.074 



2s 



1/2 



-2.07 



Total 



0.0204 



0.020 



56 A Fe 



2p 3 /2 
I/5/2 
2pi/ 2 

109/2 
24/2 
3S1/2 



-11.20 
-10.35 
-8.84 
-6.33 
-1.08 
-0.33 
Total 



0.0029 
0.0027 
0.0035 
0.0001 
0.0002 
0.0008 
0.010 



I/7/2 -5.85 
2p 3/2 -0.37 



Total 



0.0013 
0.0031 



0.004 



Table 7: Discrete contributions to the pionic decay width (units of T\) for each of the unoccupied shells, 
calculated with the NQ 7r— nucleus optical potential, as defined in the caption of Table 5. 
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The study of the mesonic decay of AA hypernuclei constitutes an obvious extension of 
this work. To improve the existing calculations [33, 37] would require a correct treatment 
of the AA pair inside of the nuclear medium. Thus, the recent work of Ref. [38], where 
short and long range correlations are taken into account, will be a good starting point for 
this end. 
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